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ABSTRACT 


Computational Fluid Dynamics (CFD) technology has been used to assess the exhaust 
plume pollutant environment of the RD-170 engine hot-firing on the FI Test Stand at Marshall 
Space Right Center. Researchers know that rocket engine hot-firing has the potential for forming 
thermal nitric oxides (NO*), as well as producing carbon monoxide (CO) when hydrocarbon fuels 
are used. Because of the complicated physics involved, however, little attempt has been made to 
predict the pollutant emissions from ground-based engine testing, except for simplified methods 
which can grossly underpredict and/or overpredict the pollutant formations in a test environment. 
The objective of this work, therefore, has been to develop a technology using CFD to describe the 
underlying pollutant emission physics from ground-based rocket engine testing. This resultant 
technology is based on a three-dimensional (3D), viscous flow, pressure-based CFD formulation, 
where wet CO and thermal NO finite-rate chemistry mechanisms are solved with a Penalty 
Function method. A nominal hot-firing of a RD-170 engine on the FI stand has been computed. 
Pertinent test stand flow physics such as the multiple-nozzle clustered engine plume interaction, 
air aspiration from base and aspirator, plume mixing with entrained air that resulted in 
contaminant dilution and afterburning, counter-afterburning due to flame bucket water-quenching, 
plume impingement on the flame bucket, and restricted multiple-plume expansion and turning 
have been captured. The predicted total emission rates compared reasonably well with those of 
the existing hydrocarbon engine hot-firing test data. 
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INTRODUCTION 


Russian-built kerosene fueled engines such as RD-170 or its likeness have been identified 
as potential candidates to fly the Single-Stage-to-Orbit Rockets. The potential of forming thermal 
nitric oxides (NO*) and of producing carbon monoxide (CO) has been a concern for ground-based 
engine testing using hydrocarbon fuels. The release of these criteria pollutants into the 
atmosphere not only contributes to acid rain and ozone depletion (NO*), but also poses as a 
potential threat (CO) to living organisms. It is therefore important to predict accurately those 
criteria pollutant emissions from engine testing for the environmental impact assessment. Several 
simplified analyses have been used in the past for first principle estimates. For example, 
thermochemical analysis using chemical equilibrium computer (CEC) code 1 can provide CO 
concentration at nozzle exit plane. However, it over-predicts CO concentration due to the 
omission of its subsequent after-burning with air. In addition, it can not predict NO formation 
unless gross assumption on the amount of air mixing can be made; A perfectly-stirred reactor 
analysis has been reported to treat the plume-air mixing 2 . Albeit simple, the condition of perfect 
mixing between the plume and entrained air is not justified. Besides, there is an inherent difficulty 
in estimating the residence time; Other methods such as steady two-dimensioal analysis of a free 
exhaust plume can be performed, However, the physics depicted is far from that occurring on a 
test stand and the evolution of the plume is a transient process. 


The rapid gains made in CFD and computer technologies have made possible the 
development of a computational methodology that can describe the pollutant emission physics 
from ground-based rocket engine testing: 3D air entrainment, 3D multiple-nozzle plume 
interaction and mixing with air, finite-rate afterburning reaction, plume impingement with flame ^ 
bucket and plume quenching through deluge water, and 3D restricted multiple plume expansion; 

In this study, a pressure-based CFD method heavily benchmarked for nozzle, plume, and 
combustion driven flows was used for this development. The hot-firing of a RD- 170 engine on 
the FI Stand was simulated. For the purpose of this study, an 11-species, 18-reaction finite-rate 
chemistry set described the after-burning. A nominal hot-firing with water-quenching 
(homogeneous two- phase formulation) was computed, along with two cases including the frozen 
chemistry and finite-rate chemistry (both without water-quenching) for comparison. The emission 
rates of the pollutants and the exhaust plume properties were computed and the effects of after- 
burning and water-quenching on the pollutant formation were compared. This resultant 
technology has potential applications in actual rocket launches and in the development of air- 
breathing engines. 


GOVERNING EQUATIONS 
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The present flow solver, Finite Difference Navier-Stokes (FDNS) code 3-4-5 is used to 
provide multi-component, steady-state and unsteady viscous flowfield solutions by solving the 
Reynolds-averaged transport equations such as mass conservation equation, Navier-Stokes 
equations, energy equation and other scalar transport equations. The general form of these 
conservation equations can be written as: 


dp<? 

dt 


+ 
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where p and q = (1, u, v, w, h, k, e and oq) stand for the fluid density and the flow primitive 
variables for the continuity, momentum, energy, turbulence model and species mass-fraction 
equations, respectively. u g stands for the grid moving speed. p e = (Pi + PtV<* represents the 
effective viscosity which is a sum of the laminar viscosity and the turbulence eddy viscosity 
divided by a turbulence modelling constant 5 , a. The source terms S q , for the governing equations 
in 3D space Xj can be written in fully conservative form as: 
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where O, Q and co, stand for the energy dissipation function, heat source and species source 
terms, respectively. G stands for the turbulence kinetic energy production rate which is written 
as: 
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(3) 


The turbulence modeling constants C ]y C 2 and C 3 are given as 1.15, 1.92 and 0.25 respectively in 

the extended k-e turbulence model 6 . The extended k-e model is superior than the standard k-e 
model 7 in that a second time scale of the production range of turbulence kinetic energy spectrum 
is added to the dissipation rate equation. This extra time scale enables the energy transfer 
mechanism of the turbulence model to respond to the mean strain more effectively. This extended 
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k-e turbulence model was rigorously benchmarked with fully developed turbulent channel and 
pipe flows, turbulent free-shear flows, flat plate turbulent boundary layer flow, turbulent flow 
over a backward-facing step, a confined turbulent swirling flow, and dump combustor flows. The 
compressibility effect on the turbulence is taken into account by using the method of Mach 
number correction 8,9 . 


SOLUTION ALGORITHM 


To solve the system of nonlinear partial differential equations, it uses finite-difference 
approximations to establish a system of linearized algebraic equations on non-staggered grid mesh 
systems. A pressure-based predictor plus multi-corrector time-marching scheme is employed so 
that flow over all-speed range can be analyzed. The time-marching scheme, total variation 
diminishing variation (TVD) discretization, and penalty function tretament of the reaction source 
terms are pertient to this work and are depicted in the following. 


Time-Marching Scheme 


The time-marching scheme is described below. For convenience, transformed equation 
(from Xj to Si system with J as the Jacobian of coordinate transformation) of Eq. (1) is written as: 


1 ap<? _ g £ . c 

j a t as, * 



(4) 


where Fj represents convection and diffusion fluxes in i-direction. First, Eq. (4) is discretized in 
finite difference form, 


-L{( M y" -(p 9 )'} = e«"'+a- 0 )R; (5) 

where superscripts n and n+1 represent old and new time levels respectively. 0 is a time marching 
control parameter and, 0 = 1.0 and 0 = 0.5 are for an implicit Euler and a time-centered time 
marching schemes, respectively. The following linearization is then incorporated. 


(p^)”' =(p«)" <«) 
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With the above approximations, the final form of the time-marching scheme can be written as: 
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dq 


W = K 


( 8 ) 


The pressure-based multi-corrector solution method is formulated using simplified 
perturbed momentum and continuity equations 3,4,5 . The simplified velocity correction equation 
can be written as: 


dpU, 

dt 


~-VP’ 


(9) 


or, in discrete form. 


-VP' (10) 

P 

p" +1 =p n +P’ (11) 


where P represents a pressure relaxation parameter (typical value of 1.0). The velocity and 
density fields in the continuity equation are then perturbed to form a correction equation. 
Neglecting higher order terms, the continuity equation can be written as. 
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Substituting Eq. (10) into (11) and letting p' = P’/RT, the following all-speed pressure correction 
equation is obtained. 
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To reduce potential oscillations in the pressure field, an upwind TVD adaptive dissipation term 
based on the density field, which is described in the next section, is added to the right hand side of 
Eq. (13). Once solution of Eq. (13) is obtained, the velocity and pressure fields are updated using 
Eqs. (10) and (11). The density field is then updated through the equation of state. The 
temperature field can also be modified by using a perturbed temperature correction equation. The 
entire corrector step is repeated 3 or 4 times such that the mass conservation condition is 
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enforced before marching to the next time level. For steady-state flow solutions, however, only 
one corrector step is used for computational efficiency. 


Second-order central differencing schemes are employed to model the diffusion fluxes and 
the source terms of the governing equations. High-order upwind schemes are used for the 
nonlinear terms, convection fluxes, to maintain solution accuracy and to enhance numerical 
stability. A third-order upwind TVD scheme 4 is employed in the present flow solver. Only the 
convection terms are modeled using the TVD flux limiters. The convection terms of the 
governing equations can be expressed by finite difference approximation as: 


— f \ i + 1/2 ft- 1/J + ^i-l/2 


where f and h represent first-order fluxes and TVD flux limiters respectively. The TVD flux 
limiters are functioned as anti-diffusion terms to recover the scheme to high-order accuracy. The 
first-order fluxes and the TVD flux limiters are given below. 


fi+in = max{o,(pt/). +I/ J<t>, + max{o,-(pt/). +[/2 }<t> i+1 


fy + 1/2 ~ 


^ |P^L 2 {^i+1/2 + ^.-1/2 +a ( ^H-l/2 ^i- 1 / 2 )}*^ > 0 

^|P^1/+l/2 {^1+1/2 ^j+3/2 + a(t/0 1+ i/ 2 — ^i+3/2 )}■ < 0 


where the minmod functions in the TVD flux limiters are written as: 

^i+l/2 = S ^8 n (^4 ) i+ 1/2 ) {o, [l^^j+1/2 1> P s *8 n (^i*i+l/2 )'^^1+1/2T 1 ]} 


The order of accuracy of this scheme is determined by the parameters a and (3. Only the second- 
order and third-order upwind schemes were used in this study. That is. 
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The compression factor, (J, is used to sharpen the contact discontinuities and slip streams for 
better wave tracking resolution. Other schemes such as second-order and third-order upwind 
schemes and a second-order central plus fourth-order dissipation scheme are also available, 
through input data selection, in the present flow solver. The option using central-difference 
scheme with artifical dissipation has been heavily benchmarked with practical rocket engine 
nozzle, plume, and combustion driven flow tests and applications such as the Space Shuttle Main 
Engine (SSME) performance and nozzle/plume flowfields (compared with JANNAF standard 
codes and hot-fire test data) 5 , SSME start-up and shut-down transients 10 , transient SSME fuel 
prebumer flow 11 . Ramjet dump combustor flows 12 , integrated combustion chamber (SSME, 
Space Transportation Main Engine, 40k) flow and heat transfer 13 , and base flow characteristics 
for a four-engine clustered nozzle configuration 14 . 


Penalty Function Treatment of the Reaction Source Terms 


For the gas-phase chemical reaction modeling, a general system of chemical reactions can 
be written in terms of its stoichiometric coefficients (Vjj and vy' ) and the i-th chemical species 
name (Mj) of the j-th reaction as 


(15) 


If we define the Arrhenius reaction rate for forward (R^) and backward (R b j), 
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the net rate of change in the molar concentration of species i due to reactions j, Ry, can be written 
as: 


08) 

and the species production rate coj (in terms of mass fraction), is calculated by summing over all 
reactions, 

<o, = S«, (19) 
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The forward reaction rate for each of the reactions is given by the modified Arrhenius law 


K f] = Aj T B/ exp(— — ) 

and the corresponding backward reaction rate is obtained using 

' K ‘, 

where Kej is the equilibrium coefficient 
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To solve the chemistry system equations numerically, an efficient penalty function 
method 15 is employed in the present study. In the penalty function method, the system of species 
equations are solved by employing a small time step size based on the assigned tolerance (0.01 
used in this study) for species mass fraction equation, that is 


A t s = min 


P ( A °0 assigned . , 

— ,i = Un 

(O, 


(23) 


To ensure species conservation and element balance conditions, the species solutions are subject 
to the following constraints: 


Z a - =L0 

(24) 

and 

0.0<a, <1.0 

(25) 

A penalty function is therefore devised to ensure the above constraints are met after every time 
marching step. This function is defined as 

h0_o| «• - a * >0.0 

a, -a* 

(26) 
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( 27 ) 


The allowable changes in species mass fractions, which are the solutions of the species continuity 
equations, are calculated by the following expression such that the second constraint, Eq. (25), is 
satisfied. 


a* +I = a* +(a* - a- ) PF (28) 

where PF = min ( PFi ), i=l, 2, ■ ■ • • n. This procedure is a crucial requirement for the 

numerical stability and accuracy of a the present model. The resulting limited changes are adjusted 
so that they are proportional to the species source terms from which the first constraint, Eq. (24), 
is satisfied. For transient reacting flow calculations, multiple steps of the penalty function 
procedures are used to reach the flow time step size based on the operator-splitting point implicit 
approach. 


BOUNDARY CONDITIONS 


The present flow solver supports various types of boundary conditions for the inlet, exit, 
symmetry, wall, cyclic/zonal boundaries and singularity lines, etc. For subsonic inlets, fixed total 
pressure or fixed mass conditions can be specified. For supersonic inlets, all flow variables are 
fixed unless transient inlet boundary conditions are required. In case of incompressible flow 
applications, only pressure field is extrapolated at the inlet boundaries. In the present application, 
the free stream boundary around the nozzles is a modified subsonic inlet with fixed total pressure 
boundary condition to allow air entrainment. This is accomplished by making two assumptions: 
1) air entrainment boundaries are sufficiently far from the nozzles such that the flow is isentropic 
and irrotational; 2) flow is allowed to entering through the air entrainment bounadries only. 
Based on assumption 1), Bernoulli's equation is used to solve for the total velocity 

q = [2(P 0 - P e )/pc]° 5 (29) 

where subscript e indicating quantities extrapolated form the interior and P G is the ambient 
stagnation pressure. Obtaining a total extrapolated velocity from the interior q e , we can compute 
a total velocity at the air entrainemnt boundary as 

q”* 1 = q e + co(q - q e ) (30) 
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where 0 ) is defined as a relaxation parameter. The flow directions are then extrapolated from the 
interior. In order not to violate assumption 2), the total velocity is set to zero if the flow 
directions indicating an outflow condition and the pressure is set to its ambient value. 

RD-170 is a regeneratively cooled, four-nozzle clustered engine which bums Kerosene 
fuel with liquid oxygen and was used to thrust Energia launch vehicles. Thermochemical analysis 1 
was performed for the thrust chamber at a nominal operating condition and the equilibrium 
products at the nozzle exit were used as the input to the propulsion system. Table 1 shows the 
computed flow properties at the chamber, throat and nozzle exit. CIT 9423 1 was used as the 
chemical formula for Kerosene fuel. 


Table 1 . RD-170 thrust chamber equilibrium flow properties 



Chamber 

Throat 

Nozzle Exit 


P, atm. 

241.96 

139.45 

0.6337 

S 

T, deg. K 

3877 

3667 

1818 

j 

i 

Mach Number 

0.000 

1.000 

3.972 

i 

1 Species Mole Fractions 




1 

CO 

0.30671 

0.30035 

0.24569 

1 

C0 2 

0.16619 

0.17879 

0.26158 

- 

1 H 

0.02245 

0.01936 

0.00016 

; 

HCO 

0.00006 

0.00004 

0.00000 

i 

H ° 2 

0.00010 

0.00006 

0.00000 


h 2 

0.07551 

0.07349 

0.09662 

i 

I H 2 o 

0.34357 

0.35801 

0.39591 

1 

I h 2 o 2 

0.00004 

0.00002 

0.00000 


0 

0.00988 

0.00741 

0.00000 


1 OH 

0.05779 

0.04820 

0.00004 


1 0 2 

0.01768 

0.01427 

0.00000 

f 


It can be seen that a significant amount of CO (24.569%) existed at the nozzle exit. This 
is the amount that could be dumped into the environemnt and can only be chemically reduced 
through afterburning. It can also be seen that there is no soot nor polycyclic aromatic 
hydrocarbon fragments produced throughout the thrust chamber, although graphite carbon and 
polycyclic aromatics were considered. This is to be expected under the nominal RD-170 
operating condition at a near stoichiometric equivalence of 1.2939. Regenerative cooling is 
another factor since there is no need for turbine exhaust dump inside the thrust chamber. 

The inlet condition to the flame trench is supersonic and fixed water mass flow rates are 
specified along the deflector wall and other wall boundaries in the flame trench. The free stream 


34 


in iii tid 




around the plume downstream of the trench outlet is frozen at see-level condition with a wind 
speed of 10 ft/sec. At exit boundaries, direct extrapolation (for supersonic outlet only), fixed 
mass or fixed pressure at a point boundary conditions can be specified. The fixed pressure 
boundary condition was used in this study for the downstream plume outlet boundary. For 
symmetry planes, which are the center planes of the computational domain, zero gradient 
boundary conditions are applied for all scalars and tangency conditions are imposed for the 
velocity vectors. Non-slip boundary conditions are employed for the momentum equations at 
solid wall boundaries. Wall-normal zero gradient pressure boundary conditions are used in the 
code. Isothermal or adiabatic wall boundary conditions can be specified separately for each wall 
segment. In the case of turbulent flow computations using high-Reynolds number turbulence 
models, wall function approaches 7,8 are employed, which is the case in the present study. 


FINITE-RATE AFTERBURNING KINETICS 


To accurately predict the contaminant concentrations of the exhaust plume, finite-rate 
chemical kinetics are included in the numerical modeling. The plume chemistry occurring in the 
flame-bucket/Test-stand flow physics includes the afterburning of CO to C0 2 , thermal NO x 
formation and decomposition, counter-afterburning effect on CO conversion due to water 
quenching and its reduction effect on NO* formation. These are described with a wet CO (12 
reactions) and a thermal NO (6 reactions) finite-rate mechanisms, as shown in Table 2. The well- 
known Zeldovich mehcanism is included in the thermal NO chemistry. For computational 
efficiency, only NO is considered in this work since other species of the No x family such as N0 2 is 
converted from NO and usually exists in trace amount. Their existence is hence included in the 
NO as “equivalent NO”. The wet CO and thermal NO mechanisms are sub-sets of a combustion 
kinetics model for complex hydrocarbon (coal derived) fuels developed for Department of 
Energy 1617-18 . Their reaction rates, in Arrhenius law form, have been validated with Jet-stirred 
combustor data, including blow-out limits, shock-tube measurements of ignition delay times 17,18 , 
and turbulent diffusion flames and flat flames data 19 . In addition, the thermal NO x reactions have 
also been benchmarked with an industrial burner data 20 . 

The formation of thermal NO is significantly influenced by flame temperature than other 
types of NO, e.g., the fuel NO resulting from nitrogen compounds contained in the fuel. Its 
strong dependence on temperature results from both the temperature dependence of the forward 
rate constant of reaction O + N 2 = N + NO and the sensitivity of O atom equilibrium 
concentration to temperature. Production of thermal NO is generally negligible at low 
temperatures. It is therefore expected that most of the thermal NO will be formed in the flame 
front, .i.e., the plume mixing layer near the exit plane of the nozzles. 
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Table 2. Afterburning chemical kinetics 


I" Reaction 1 

A 

B 

E/R 





H 2 + 0 2 = OH + OH 

1.7000E13 

0 

2.4070E4 

OH + H 2 = H 2 0 + H 

2.1900E13 

0 

2.5900E3 

OH + OH + 0 + H 2 0 

6.0230E12 

0 

5.5000E2 

0 + H 2 = H + OH 

1.8000E10 

1.0 

4.4800E3 

H + O 2 - O + OH 

1.2200E17 

-0.91 

8.3690E3 

M + 0 + H = OH + M 

1.0000E16 • 

0 

0 

m + o + o = o 2 + m 

2.5500E18 

-1.0 

5.9390E4 

m + h + h = h 2 + m 

5.0000E15 

0 

0 

M + H + OH = H 2 0 + M 

8.4000E21 

-2.0 

0 

CO + OH = H + C0 2 

4.0000E12 

0 

4.0300E3 

CO + o 2 = co 2 + 0 

3.0000E12 

0 

2.5000E4 

1 CO + 0 + M = CO 2 + M 

6.0000E13 

0 

0 





0 + N 2 = N + NO 

1.3600E14 

0 

3.7750E4 

N 2 + 0 2 + NO + NO 

9.1000E24 

-2.5 

6.4600E4 

NO + 0 = 0 2 + N 

1.5500E9 

1.0 

1.9450E4 

M + NO = 0 + N + M 

2.2700E17 

-0.5 

7.4900E4 

N + OH = NO + H 

4.0000E13 

0 

0 

C0 2 + N = CO + NO 

2.0000E11 

-0.5 

4.0000E3 






* M stands for third-body collision partner 


COMPUTATIONAL GRID GENERATION 


The grid generation for the nozzles and test stand was performed using an EZSURF 
code 21 . It was used to interactively create the edge curves of the nozzle exits, aspirator, flame 
deflector and multi-zone block edges. The initial surfaces were then generated using transfmite 
interpolation in EZSURF. The flame deflector and nozzle exit surfaces were then elliptically 
smoothed. Further work was done on the nozzle exit surfaces using Bezier curve and local 
redistribution in EZSURF. The volume grid for the first block (Zone 1) was created using two 
linear stackings; one from the top of the block to the nozzle exit plane and then another form the 
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nozzle exit plane to the bottom of the aspirator. The flame deflector block (Zone 2) and 
subsequent external ambience (Zone 3) volume grids were created using transflnite interpolation. 

The FI Test Stand, standing 230 feet tall with a flame bucket (deflector) attached to the 
aspirator, was used to test FI engines with which Saturn launch vehicles were propelled. Not 
only does the flame bucket quenches the rocket exhaust plume with deluge water, but also turns 
the vertical flowing exhaust plume to that of a horizontal direction, after which the plume expands 
and dissipates into the atmosphere. Fig. 1 shows the computational domian for the F-l Stand. 
The RD-170 engine is mounted vertically, firing down into the flame bucket. Due to symmetry, 
only half of the domain was actually computed. The four RD-170 nozzles (mounted beneath the 
platform that is not modeled) and the aspirator are described by Zone 1 which contains 63,360 
grid points (72x40x22). The aspirator itself is mounted on top of the flame bucket such that the 
air entrainment can be promoted and the plumes are centered while impinging at a predetermined 
area in the flame bucket (approximately 45-degree elbow at the bottom). Ambient air is allowed 
to be entrained through the top and four side boundaries of Zone 1. The flame bucket is modeled 
by Zone 2 which composes of 72,000 grid points (72x40x25). The plumes are then quenched 
through water deluge injecting from all four walls inside the bucket. The water injection pattern is 
designed as such that most of the water injects at the plume impingement area. 

After passing through the flame bucket (Zone 2), the quenched plumes expand into the 
vast surrounding atmosphere (Zone 3) which is described by 156,975 grid points (91x69x25). 
The relative sizes and loactions of all three zone are shown in Fig. 2. Total number of grid points 
used in this study was 292,335. 


RESULTS AND DISCUSSION 


A frozen chemistry analysis was performed at first for a small period of elapsed time in 
order to establish an initial plume in the system. This not only prepared for the initial flowfields 
for the computation of subsequent parallel finite-rate chemistry and finite-rate chemistry with 
water-quenching cases, but also served as an excellent check of mass conservation using carbon 
balance. This procedure is allowed since the total pollutant emission rates during a nominal 
steady RD-170 hot-firing are desired and not those of an actual start-up sequence. The goals are 
therefore to compute the growth rate of NO and the disappearance rate of CO until they achieve 
asymptotic states. 

Fig. 3 shows the computed velocity vectors colored by Mach number contours to 
represent the entrained air surrounding the bulk rocket plumes. Entrained from surroundings 
close to the open platform, the air accelerated and mixed with the plume bounadry layer and 
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entered the opening of the aspirator, where it continuously mixed with the advancing plume 
bounadry layer. In addition to cooling and dilution of the contaminant, the entrained air also 
provides the source of reactants for afterburning and thermal NO formation. Due to a backward- 
facing step geometry between the aspirator and the flame bucket, flow recirculation patterns 
formed underneath the aspirator. 

Fig. 4 shows the computed temperature contours for the water-quenching case. Water jet 
vectors are shown for the top and bottom walls. The water jets at the back wall are mostly 
blocked out by the center temperature contours, whereas the water jets from the front wall are not 
shown for clarity. Most of the thermal NO is formed near the aspirator level where it has the 
most mixing and the hottest temperature. This is in agreement with the characteristics of the 
Zeldovich mechanism. It can also be seen that the plumes impinge on the 45-degree bend section 
of the flame bucket where it encounters the most water-quenching, that agrees with the original 
design. The quenched plumes then turn and partially hit the outer wall, where it moves 
horizontaly out, following the direction of the flame bucket. 

Fig. 5 shows the comparisons of averaged mass fraction in the system for species CO, 
C02, NO and OH, with respect to elapsed time. The averaged mass fraction for CO in the add 
water-quenching case is more than that of the purely finite-rate chemistry case due to the counter 
after-burning effect of deluge water, and vice versa for that of C02. The higher level of OH in 
the finite-rate case shows a higher degree of after-burning reaction, due to higher overall plume 
temperature without water-quenching. Total NO production drops significantly in the add water- 
quenching case. The concentraion of NO in the add water-quenching case is almost two orders of 
magnitude lower than that in the finite-rate chemistry only case. This is not surprising since the 
extent of thermal NO formation depends heavily on the local temperature. The effect of water 
deluge on the formation of NO is the reduction in peak temperature caused by the heat capacity of 
water. 


The computed growth of plume volumes is shown in Fig. 6. It can be seen that the 
growth rates of the plumes have reached their approximate asymptotic states. Obviously, the 
plume energy of the finite-rate chemistry is higher than that of the add water-quenching case. 
Correspondingly, the plume volume growth rate of the finite-rate case is larger than that of the 
add water-quenching case. It is anticipated that the characteristics of computed growth of plume 
energies would be similar to those of plume volumes. These CFD results ultimately serve as the 
basis (input) for the subsequent meteorological cloud dispersion calculation, where the plume 
volume growth rate helps determining the eventual plume size and the plume energy growth rate 
helps determining the magnitude the plume buoyancy force. 

A comparison of the calculated criteria pollutant total emission rates for RD-170 with 
those measured for other Kerosene-fueled engines 22 is shown in Table 3. Since thrust levels are 
quite different, the measured emission rates were extrapolated to a RD-170 level by thrust ratios. 
Although the operations of the engines and test stands are vastly different, and there is a question 
on whether the point sampling technique used in the measurement representative of the whole 
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plume. The agreement in terms of order-of-magnitude for these engines is reasonable and 
encouraging. Among RD-170 CFD computations, the emission rate of CO is maximum and that 
of NO is zero for the frozen chemistry case, since afterburning reactions are not turned on, 
whereas the effect of water deluge has reduced the CO conversion rate from 83% of the finite- 
rate chemistry case to that of 67%. This is of interest since the effect of water deluge also has 
reduced the NO x production by 95%. 


Table 3. Total emission rates 


Engine 

Thrust* Ibf 

NO x , lb/sec 

CO, lb/sec 

MA5B/hot-firing 

370,000 

5.4/25.9’ 

133/641’ 

MA3S/hot-firing 

165,000 

2.7/29.1* 

210/2,266’ 

MA3B/hot-firing 

60,000 

1.5/45.6’ 

138/4,111’ 

RS27/hot-firing 

205,000 

1.2/10.6’ 

94/820’ 

RD-170/CFD 




Frozen 

1,777,000 

- 

1,382 

Finite-Rate 

1,777,000 

8.0 

232 

Add Water 

1,777,000 

0.4 

463 


* based on extrapolation of measured emission rate to a RD-170 by thrust ratio 


CONCLUSION 


A 3D viscous flow, pressure based CFD technology has been developed to predict the 
contaminant emissions from ground-based rocket engine RD-170 propulsion testing. Pertinent 
test stand flow physics such as the multiple-nozzle clustered engine plume interaction, aspiration 
from base and aspirator, plume mixing with entrained air that resulted in contaminant dilution and 
afterburning, counter-afterburning due to flame bucket water-quenching, plume impingement on 
the flame bucket, and restricted multiple-plume expansion and turning, have been predicted. The 
predicted criteria pollutant total emission rates agreed reasonably well with those of the existing 
hydrocarbon engine hot-firing test data. This resultant technology has potential applications in 
actual rocket launches and in the development of air-breathing engines. 
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Fig. 1 Computational domain for the RD-170 nozzles and F-l Test Stand. 
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Fig. 2 Computational domain for all three zones. 
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Fig. 3 Velocity vectors colored by Mach number for the air entrainment 



CFD Simulation of Kerosene Engine Test at MSFC 

(Water-quenching in Flame Deflector) 



Fig. 4 Computed temperature contours for the water-quenching case. 
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Fig. 6 Comparison of plume volume 






